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Abstract: The annihilation or decay of Dark Matter (DM) particles could affect the 
thermal history of the universe and leave an observable signature in Cosmic Microwave 
Background (CMB) anisotropics. We update constraints on the annihilation rate of DM 
particles in the smooth cosmological background, using WMAP7 and recent small-scale 
CMB data. With a systematic analysis based on the Press-Schechter formalism, we also 
show that DM annihilation in halos at small redshift may explain entirely the reionization 
patterns observed in the CMB, under reasonable assumptions concerning the concentration 
and formation redshift of halos. We find that a mixed reionization model based on DM 
annihilation in halos as well as star formation at a redshift z 2^ 6.5 could simultaneously 
account for CMB observations and satisfy constraints inferred from the Gunn-Peterson 
effect. However, these models tend to reheat the inter-galactic medium (IGM) well above 
observational bounds: by including a realistic prior on the IGM temperature at low redshift, 
we find stronger cosmological bounds on the annihilation cross-section than with the CMB 
alone. 
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1. Introduction 

In the minimal ACDM model, the CMB has very little to say about Dark Matter (DM), 
apart from a measurement of the relic abundance parameter Q^j^yih^ . However, there is a 
chance that DM could leave another signature in the CMB. In the case of annihilating DM, 
if the ratio of the annihilation cross section over the mass is not too small, annihilation 
products could contribute to the ionization of the thermal bath, and affect the history 
of recombination and reionization. This has already been discussed in detail in several 
references including [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. In the case of decaying dark matter, 
similar effects could take place if the particle lifetime is not too large [13, 14, 15, 16, 17]. 
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In a detailed study of these mechanisms, assumptions concerning the nature of dark 
matter are of course crucial. Different DM particles may annihilate or decay in different 
channels (into hadrons, leptons, gauge bosons, etc.). The produced particles can themselves 
decay in several steps, until only stable particles like electrons, photons and neutrinos re- 
main. While the energy contained in neutrinos is lost for the cosmic plasma, other decay 
products can contribute to the ionization and heating of this plasma. The authors of [7] 
computed /(^), the redshift-dependent fraction of the energy produced by annihilations 
that contributes to the ionization and heating of the plasma, for several WIMP models. 
Similar calculations could be carried on any type of annihilating or decaying DM particles. 
In a cosmological analysis like that presented in this paper, we don't need to refer specifi- 
cally to a given DM particle, and we may incorporate all model-dependent particle physics 
assumptions in the unknown fraction f{z). 

In this work, we revisit the impact on the CMB of annihilating DM. We will not 
introduce any new physical ingredients, neither from the point of view of particle physics, 
nor from that of structure formation models. However, we will present some new generic 
parametrizations of the relevant phenomena, in order to perform a systematic comparison 
of various models to recent CMB data sets. This will allow us to tighten some of the bounds 
presented previously in the literature. 

In section 2, we summarize the impact of DM annihilation on recombination, and 
explain how we took it into account by modifying the public Boltzmann code CLASS. This 
code can simulate recombination with either of the two public codes RECFAST and HYREC. 
We double-checked our results by modifying the two algorithms. We show that they give 
the same results, but for non-trivial models like those assumed in section 4, only HYREC 
remains numerically stable. 

In section 3, we focus on the effects of annihilation in the smooth DM background 
distribution. We find similar but slightly stronger bounds than in recent studies, thanks 
to our updated CMB data set. With a generic parametrization of the redshift-dependence 
of the function accounting for the fraction of energy released to the gas, we confirm that 
current CMB data is not sensitive to this dependence. 

In section 4, we consider the additional effect of enhanced DM annihilation in halos at 
small redshifts. This effect has been previously discussed in several references including [3, 
4, 5, 8, 11, 12]. Some of these works suggest that it could account for a significant fraction (if 
not the totality) of the reionization of the universe at low redshift. We derive an approximate 
but rather generic parametrization of this effect, and carry on the first systematic parameter 
inference using current CMB data in a model with reionization from annihilation. We also 
confront to the data a mixed model, with reionization explained both by DM annihilation 
and star formation, and discuss the relevance of this model for explaining simultaneously 
CMB data and Gunn-Peterson bounds. 

Our general conclusions and future directions of research are outlined in section 5. 
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2. Equations of recombination with Dark Matter annihilation 



2.1 Standard recombination 

Before discussing the impact of DM annihilation (or alternatively dark matter decay) on 
the recombination history, we first briefly review the standard recombination model. We 
start by describing the simple three-level atom model of Peebles [18, 19] and then discuss 
the subsequent improvements of this model. 

In what follows we shall assume that helium has entirely recombined (which is indeed 
the case for all redshifts of interest) and only deal with hydrogen recombination. We denote 
He the number density of free electrons, uh the total number density of hydrogen (in ionized 
and atomic forms), — ^ the free electron fraction and Tm (resp. T^) the matter (resp. 
photon) temperature. 

2.1.1 The effective three-level atom (TLA) model 

It is well known that direct recombinations to the ground state are highly ineflicient: if 
a hydrogen atom forms directly in its ground state, it emits a photon which is going to 
immediately ionize an other atom, leaving the overall free electron fraction unchanged. 
The basic idea of Peebles' "case B recombination" is that efficient recombination only takes 
place when the electron gets first captured into an excited state n > 2, from which it 
cascades down to n = 2. The newly formed atom may then eventually reach the ground 
state, either by emitting a Lyman-a photon from 2j), or by the 25 Is two-photon process. 
The Lyman-a line being very optically thick, the net rate of 2p Is transitions is, to a 
first approximation, the rate at which Lyman-a photons redshift across the resonance. At 
early times {z > 900), the net rate of transitions from the n — 2 state to the ground state is 
much smaller than the rate at which excited atoms are photoionized by CMB photons, and 
the slow 2^1 transitions constitute the bottleneck of the recombination process. At late 
times [z < 900), the intensity of the radiation field drops, and atoms that do recombine 
to an excited state almost certainly reach the ground state; during this period the new 
bottleneck is the rate at which free electrons and protons can encounter each other and 
recombine. 

To correctly describe recombination, accounting for the effects mentioned above, Pee- 
bles introduced the pre-factor C defined as 

where Ah = 8.22458 s"-*^ is the decay rate of the 2s level, Kh = s7tH{z) accounts for the 
cosmological redshifting of Lyman-a photons^, and i3h is the effective photoionization rate 
from n = 2 (per atom in the 25 state). C represents the probability for an electron in the 
n = 2 state to get to the ground state before being ionized. The evolution equation for the 
free-electron fraction is then given by 



dx^ 1 
~dz ^ {l + z)H 



(2.2) 



^[KHTiHi^ — Xe)] ^ is the rate of escape of Lyman-a per atom in the 2s state 
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where an is the case-B recombination coefficient and z/q, is the Lyman-a frequency. Because 
the effective recombination rate per free electron CanXeriH is always much smaller than 
the Hubble rate (due to the two bottlenecks mentioned above), primordial recombination 
proceeds much slower than in Saha equilibrium. 

In addition, the matter temperature is determined from the Compton evolution equa- 
tion: 

dTu _ 1 SardrT^ Xe ^^^^ - T ) + 
dz (1 + z)H SrrieC 1 + fne + ^ I + z 

where we have defined the dimensionless parameter 

_ SaTdrT^ Xe 
^ ~ 3HmeC 1 + fne + ' 

where ay is the Thomson cross-section, the radiation constant, rUe the electron mass, 
c the speed of light and /hc the fraction of helium by number of nuclei. When 7^1, 
the matter temperature is locked to the radiation temperature by Compton heating, Tm ^ 
oc (1 + z); it decays adiabatically as Tm oc (1 + z)^ when 7^1. 

2.1.2 Beyond the TLA model 

With the prospect of upcoming high-precision data from the Planck satellite, several groups 
have revisited the simple TLA model presented above and introduced important corrections. 
Here we use the codes RECFAST [20] and HYREC [21] which implement these corrections, 
approximately for the former and exactly for the latter. The corrections are of two types: 

• Highly excited states of hydrogen are not in Boltzmann equilibrium with each other, 
and one must account for all bound-bound and bound-free transitions involving them, in- 
cluding stimulated transitions. The original RECFAST code accounted for these transitions 
approximately by multiplying the case-B recombination coefficient (and effective photoion- 
ization rate) by a "fudge factor" F = 1.14, fitted to reproduce multilevel computations 
[20]. The original computation of Seager et al assumed that angular momentum substates 
were in statistical equilibrium. This approximation, however, was shown not to be accurate 
enough [22, 23]. The latest version of RECFAST attempts to account for these more detailed 
high-n computations approximately by using a new fudge factor F = 1.105. 

It turns out that the effect of highly-excited states can be exactly and efficiently ac- 
counted for by generalizing the case-B coefficient to a non-zero CMB temperature [24, 25]. 
The code HYREC is using precomputed effective recombination coefficients in an effective 
few-level atom model, at virtually no speed cost compared to the TLA model, and with the 
advantage of being exact. 

• Being the recombination bottleneck at early times, the Lyman-a escape and two- 
photon 25 Is decays need to be modeled very precisely. Several radiative transfer effects 
were shown to be important for high-accuracy predictions of CMB anisotropics (see for 
example Refs. [26, 27, 28, 29, 23] and references therein). Detailed codes such as HYREC 
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(2.3) 
(2.4) 



and COSMOREC [23] account for all important radiative transfer effects exactly, by evolving 
the radiation field numerically. This part of the calculation is heavier computationally, but 
efficient implementations render the runtime for the recombination calculation comparable 
with the runtime of the Boltzmann code itself. RECFAST accounts for radiative transfer 
effects by adding a correction function to the recombination rate Xelcorr, fitted to reproduce 
the detailed codes for cosmologies close to the current best-fit value. 

2.2 Parametrization of Dark Matter annihilation 

We wish to express the rate at which the energy released by DM annihilations is injected in 
the thermal bath. In the next subsection, we will summarize how this energy is used and 
in which proportions. 

We can write the energy injected into the plasma per unit of volume and time as the 
product of the number of DM particle pairs ripairsi the annihilation probability per unit of 
time Pann, the released energy per annihilation E'ann, and the redshift-dependent fraction 
of released energy f{z) absorbed by the gas 

av at z 

= pynU{l + zfm^. (2.5) 

In Eq. (2.5), a is the annihilation cross-section, v is the relative velocity of DM particles, 
{av) is the average of axv over the velocity distribution, moM the mass of the DM particle, 
pc = the critical density of the universe today (with Hq the Hubble constant today), 
and f^DM the Dark Matter abundance today relative to the critical density. In the case 
where DM consists of Dirac Fermions, there should be an extra factor 1/2 in the last two 
equalities (since only half of the pairs are made of one particle and one anti-particle). If 
this is the case, we can decide to absorb this factor in a redefinition of f{z). Then, for a 
given cosmological evolution, all the model-dependent part of the energy injection rate can 
be parametrized by the following function of redshift, 

P.nn{z) = f{z)^. (2.6) 

The authors of [7] computed f{z) for several WIMP models. They found that f{z) is a 
smoothly decreasing function, lying in the range from 0.2 to 0.9 at redshift 2500 (depending 
on the WIMP mass and dominant annihilation channel), and decreasing by a factor 2 to 
5 at small redshift. Similar calculations could be carried for any type of annihilating or 
decaying DM particles. 

2.3 Effects of Dark Matter annihilation on the thermal history of the universe 

The energy injected by DM annihilation has three effects: ionizing the plasma, exciting 
hydrogen atoms , and heating the plasma [30, 6]. A fraction of the atoms excited by the 
second mechanism will be subsequently ionized by CMB photons. Hence, the first two 
effects (illustrated in Figure 1) have a direct impact on the free electron fraction, and the 
last one on the matter temperature. 
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In this section, we will use the on-the-spot 
approximation, which supposes that all inter- 
actions between the decay products of DM an- 
nihilation and the plasma take place locally, 
on a time scale negligible with respect to the 
expansion time scale. It can be shown that 
this approximation is very good for z ^ 100 
(see section 4.2 and figure 7). In section 3, we 
will only be interested in annihilation effects 
at such high redshifts, and we will stick to this 
approximation. In section 4, for the study of 
low-redshift reionization, it will be necessary 
to go beyond. 

Under such an approximation, the three 
effects of DM annihilation can be described 
with the following equations: 



+ e- 




Figure 1: Illustration of the impact of DM 
annihilation on Peeble's "case B recombina- 
tion" model. 



Ionization of the plasma. In the on-the-spot approximation, if Xi(^) denotes the fraction 
of the injected energy going into ionization, and Ei the average ionization energy per baryon, 
the number of direct ionizations per interval of redshift dz reads 



dE 



{1 + z)H{z)nH{z)Ei dVdt 



(2.7) 



DM 



where we used dz/dt — —(1 + z)H. Shull and Van Steenberg found that for a neutral gas, 
approximately 1/3 of the energy goes into ionization [31]. For an ionized gas, none of the 
energy can be used for ionization. Thus, for a partially ionized gas, Chen and Kamionkowski 
proposed to approximate Xi by (1 — a;e)/3 [32]. The fact that Xi (1 — ^e) makes physical 
sense, since the ionization rate must be proportional to the abundance of neutral hydrogen. 



Excitation of hydrogen. The rate of collisional excitation of hydrogen (ls^2p and 
ls^2s transitions, etc.) due to DM annihilation is similar to that of direct ionization, with 
Ei replaced by the Lyman-a energy E^ and Xi(^) by the fraction Xa{z) of the injected 
energy going into excitations. Once a given atom is in the n — 2 state, it has a probability 
{1 — C{z)) to be ionized by CMB photons. Thus, the net ionization rate per redshift interval 
dz due to collisional excitations by DM annihilation products reads 



dE 



{l-C{z))Xa{z) 

{1 + z)H{z)nH{z)Ea dVdt 



(2. 



DM 



Chen and Kamionkowski showed that in first approximation one may assume Xi — Xn — 
(1 — Xe)/3. Note that this process is subdominant with respect to the direct ionization of 
the plasma. 
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Heating of the plasma. Finally, DM matter annihilation heats the plasma at a rate (per 
unit of time) 



dTu 



dt 



2Xh 



dE 



(2.9) 



DM 



DM 3ki,nH{^ + fne + Xe) dV dt 
with Xh = 1 — Xi ~ Xq; = (1 + 2xe)/3 the remaining fraction of the total injected energy. 
2.4 Recombination equations with DM annihilation 

We can now write the modifications needed for each of the two recombination codes RECFAST 
and HYREC, both implemented in the Boltzmann code CLASS^ [33, 34] used throughout this 
work. The point of using two different codes is to compare the results and check that our 
approach for including annihilation effects is robust and consistent. In addition, we will see 
that in some of the cases discussed below, the second code is more stable numerically and 
allows to explore more general models. 

Given equations (2.5 - 2.9), implementing DM annihilation in the two codes only 
requires to add two new terms proportional to Pann 

(z) in the basic equations for hydrogen 

recombination and matter temperature: 







dz 


dz 


dTu 


_ dTu 


dz 


dz 



H{z) 



1 - Xe{z) 



C{z) 



St 



H{z) 



2 l + 2xe{z) 



1 



Sfcft "iriHiz) l + fHe + Xe{z 



(2.10) 

(2.11) 



where the subscript "st" stands for the standard rates, given by eqs (2.2,2.3) for case B 
recombination. These equations neglect the possibility that a fraction of the energy released 
by DM annihilation would serve for helium ionization. As in ref. [9], we checked that such 
a refinement would have a negligible impact on the CMB spectra. 

In HYREC, it is also necessary to write separately the impact of DM annihilation in 
the equations accounting for approximation schemes, like the steady-state approximation 
for the matter temperature at early times. In appendix A, we write explicitly our modified 
HYREC equations. 



In section 4, we will introduce extra modifications allowing to account for DM annihi- 
lation at small redshift beyond the on-the-spot approximation. Our modification to CLASS, 
including those in RECFAST and HYREC, will be part of the next public distribution 1.5 of 
the code. 



3. Dark Matter annihilation before structure formation and reionization 

In this section, we wish to better understand the impact of DM annihilation on the CMB 
at relatively high redshift, i.e. roughly for z ^ 100. At lower redshift, enhanced DM 
annihilation in non-linear structures might be responsible for additional effects that we will 

^http : / / class-code . net 
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study separately in the next section. Since the two regimes have a rather different impact on 
the CMB spectra, it is legitimate to split the discussion in this way. DM annihilation effects 
on the CMB at high redshift have been thoroughly investigated by Galli et al. [6, 9, 10]. 
In this section, we will only update previous results, before exploring new models including 
halo effects in the next section. 

For simplicity, we first assume in subsections 3.1, 3.2, 3.3 that the annihilation pa- 
rameter j)ann IS independent of redshift, as in [6, 9]. We will relax this assumption in 
subsection 3.4. 

3.1 Annihilation effects on and Tm 

In figure 2, we show the evolution of Xe(z) and Tm{z) computed with either RECFAST or 
HYREC for four values of the annihilation parameter. We tested RECFAST and HYREC 




10° 10^ 10^ 10^ 10^ 10^ 10° 10^ 10^ 10^ 10^ 10^ 



1+z 1+z 

Figure 2: Free electron fraction and matter temperature as a function of the redshift with, from 
bottom to top, Pann = 0, 10~^, 5- 10~^ or 10~^ m^s~^kg~^. For each value of Pann, we used either 
RECFAST or HYREC, and two different options for each of the two codes; the four results agree to 
better than a few percent, and the difference would be indistinguishable on the plots. 

in two modes: for RECFAST, with or without taking into account the hydrogen physics ef- 
fects described in [35] (using the switch Hswitch), and for HYREC, using the mode RECFAST 
(mimicking a simplified version of RECFASt) and FULL (including a state-of-the art descrip- 
tion of an effective multi-level hydrogen atom as well as radiative transfer near the Lyman 
lines). The FULL mode uses interpolation tables requiring Tm < 7^. This is the case at all 
times provided that the annihilation parameter does not exceed j^ann ^ 3- 10~^ m^s~^kg~^. 
In order to test HYREC/FULL above this value, we removed the condition Tm < from the 
code, letting it extrapolate from the table. For all used values of Pann, never exceeds 
by a large fraction and the extrapolation is therefore accurate. 

In the results presented in figure 2, we assumed a ACDM model without reionization. 
The first two small steps seen on the electron fraction curve correspond to the two helium 
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recombinations, and bring the ratio Xg = n^jna down to one. The third and biggest step 
accounts for hydrogen recombination. As expected, the energy injected by DM annihilation 
inhibits recombination, and the free electron fraction freezes out at a larger value. More- 
over, the matter temperature decreases more slowly after photon decoupling due to energy 
injection in the gas resulting from DM annihilation. 

For each value of pann, the difference between the four algorithms is extremely small. We 
checked that the shifts induced in the CMB power spectra are well below the sensitivity level 
of current CMB data sets, and lead to the same observational bounds on Pann- This means 
that the four approaches can be used indifferently in the rest of this analysis. Whenever we 
could, we sticked to RECFAST with Hswitch on, in order to speed up the computation. We 
will mention below that for some models, we had to use instead HYREC with the RECFAST 
or FULL mode, found to be the more stable numerically. In these cases, the increase in 
computing time in the full parameter extraction process was less than a factor of two. 

3.2 Effects on the CMB Power spectrum 

We could expect the effect of DM annihilation to be degenerate with that of reioniza- 
tion, since both mechanisms increase the ionization fraction after photon decoupling, and 
therefore the optical depth to last scattering t{z^q^. Indeed, a high ionization fraction at 
z < Zdec implies that more photons interact along the line of sight, which tends to damp 
temperature and polarization anisotropics on sub-Hubble scale, and to regenerate extra 
polarization around the Hubble scale at the time of re-scattering. 

In figure 3, we compare the effect of varying pann with that of changing the redshift of 
reionization, under the usual simplifying assumption of a single reionization step, such that 
Xe{z) follows a hyperbolic tangent centered on Zreio- The two effects turn out to be rather 
different for reasons that are easy to understand. 

First, the annihilation effect is already present around z = z^ec? and results in a small 
delay in the decoupling time (defined as the maximum of the visibility function —r^e~^). 
Hence, the sound horizon at decoupling has the time to grow, while the diffusion damping 
scale has sufficient time to reach larger scales. The increased sound horizon results in peaks 
visible under larger angles or smaller /'s: this shifting of the peak explains the oscillatory 
patterns clearly visible in figure 3. The increased diffusion damping scale enhances Silk 
damping at large /'s, leading to the negative high-/ slope in figure 3. 

Second, DM annihilation increases the ionization fraction and the optical depth at 
all redshifts in the range < z < Zdec- This means that some power is removed from 
the temperature and polarization spectrum on all scales, with a maximum suppression 
for / > 200, corresponding to modes being always inside the Hubble radius in the range 
< z < Zdec- In the temperature spectrum, multipoles with / < 200 are less and less 
affected when / decreases. In the polarization spectrum, the rescattering of the photons 
generate extra polarization on all scales in the range 2 < / < 200 corresponding to the 
variation of the Hubble scale between decoupling and today. 

In contrast, reionization enhances Xe{z) only at small redshift, z < 10 in our example. 
It does not affect recombination and does not shift the peaks: the reionization curve in 
figure 3 has no oscillatory patterns. Power is maximally suppressed for all scales being 
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10' 10^ 10^ 




10^ 10^ 10^ 

multipole I 



Figure 3: Temperature and polarization power spectra for several models with DM annihilation or 
reionization, rescaled by a reference model with none of them. The curves with oscillatory patterns 
correspond to different values of Pann (expressed in the key in m^s~^kg~^) and no reionization. The 
last curve was obtained with Pann — and with reionization at z^eio — 11. 

inside the Hubble radius at reionization, i.e. all multipoles / > 20. The regeneration of 
power in the polarization spectrum is limited to / < 20 for the same reason (but is very 
strong, since reionization enhances Xe{z) much more than DM annihilation). 

Hence, DM annihilation effects are clearly not degenerate with reionization effects. In 
order to check that the impact of pann cannot be mimicked by other parameters in the 
ACDM model, we should however run a parameter extraction code and marginalize the 
posterior distribution of pann over other cosmological parameters. 

3.3 Analysis with WMAP and SPT data 

We compared to observations a model described by the six free parameters of the vanilla 
ACDM model, the annihilation parameter Pann, and the effective neutrino number A^eff? 
accounting e.g. for extra relativistic degrees of freedom. Since the South Pole Telescope 
(SPT) collaboration reported an intriguingly high best-fit value of A'eff [36], we wish to check 
whether pann and A'eff are correlated in some way, such that the effect of one parameter 
could be confused with that of the other. A priori, this is not impossible, because both 
parameters impact the amplitude of the high-/ damping tail of the temperature spectrum, 
relatively to the amplitude of the first acoustic peaks. 

We compared this model to WMAP 7- year data [37] and SPT data [36], using the 
code MONTE PYTHON [38], based on Monte Carlo Markhov Chains and on the MetropoHs- 
Hastings algorithm (Hke COSMOMC [39], but MONTE PYTHON is interfaced with CLASS 
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instead of CAMB [40], written in PYTHON, and has extra functionalities; this software wih 
soon be released publicly). On top of the cosmological parameters, we vary three nuisance 
parameters related to the foreground contamination of the SPT data and constrained by 
gaussian priors, following strictly the recommendations and the software released by the 
SPT collaboration. All results on cosmological parameters are marginalized over these three 
nuisance parameters. We took flat priors on all parameters and just imposed pann > 0. 

Our results, summarized in the flrst column of Table 1 and in the triangle plot of 
figure 4, are in excellent agreement with those of the SPT collaboration for the first seven 
parameters (last column of Table 3 in ref. [36]). For DM annihilation, we obtain a bound 

Pann < 0.44 X 10"VVs/kg (WMAP7 + SPT, 95%C.L). (3.1) 

We observe no correlation between pann and any other parameter in the analysis (in partic- 
ular, we checked that there is no correlation at all with A^eff)- The marginalized posterior 
probability for 2?ann is displayed on the left plot of figure 4, and shows no evidence for DM 
annihilation in current data. Our bound is stronger than the most recent one, presented in 
ref. [9], 

Ps^nn < 2.09xlO-2^cmVs/[GeV/c2] = l.lTxlQ-VVs/kg (WMAP7 + ACT, 95%C.L), 

(3.2) 

due to the inclusion of the SPT dataset. It is also stronger than that from ref. [11]. 
We refer the reader to ref. [9] for a discussion of derived limits on the DM annihilation 
cross-section, given various ansatz for the mass and released energy fraction /. 

3.4 Redshift dependent annihilation parameter 

In any realistic model, the fraction of energy absorbed by the overall gas is a function of the 
redshift / = /(^), as shown in Figure 4 of Slatyer et al. [7] for several examples. The shape 
of f{z) depends on the DM annihilation channel(s). The impact of the redshift-dependence 
of f{z) on current /future CMB constraints on DM annihilation has been questioned with 
different methods in various papers [11, 9, 10]. Here we will check this issue with yet another 
approach, and confirm the results of other references showing that taking this dependence 
into account is of very little relevance. 

All examples for f(z) shown in [7] have strong similarities: f{z) is always a smooth 
step-like function, with plateaux at z > 2500 and z < 30. In view of performing a model- 
independent comparison to the data, it is tempting to approximate f{z) with a family of 
simple analytic functions, capturing the essential behavior of f{z) in all cases. The CMB is 
marginally affected by the behavior of f{z) at low z even in the case of a constant Pann, and 
even more if f{z) decreases; moreover, the effect of DM annihilation at low z is superseded 
by that of reionization. Hence, a given approximation scheme doesn't need to be accurate 
at low z, but should capture the essential behavior for z > 100. Figure 4 in [7] suggests 
that f{z) could be chosen to be constant at z > 2500, to decrease like a parabola in log- log 
space for 30 < z < 2500, and to remain again constant at z < 30, as displayed in figure 5. 
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Table 1: Mean and edges of the 68% Minimum Credible Interval (MCI) for the cosmological 
parameters of the five models that we compared to WMAP7 and SPT data. We don't show results 
for the nuisance parameters associated to SPT data, that have been marginalized over. The second 
model differs from the first one by the inclusion of a 2:-dependent annihilation function parametrized 
by a. All parameters have been assigned top-hat priors, and never reach prior edges except Pann 
(limited to positive values), a (limited to the range —0.2 < a < 0) and (on which we imposed a 
prior 20 < < 30 only in the last column). For pann (and /h in the last column), we indicate the 
95% Confidence Level (C.L.) upper bound. 

The work of [10] analyzed the amount of information that one can extract from current 
and future CMB data on f{z) or Pann(^) (we recall that these two functions are simply re- 
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Figure 4: One and two-dimensional marginalized posterior probabilities for the free parameters 
of a ACDM model with a free effective neutrino number A/eff and time-independent annihilation 
parameter Pann, compared to WMAP7 + SPT data. For the two-dimensional posterior, we show the 
contours corresponding to the 68.3%, 95.4% and 99.7% credible regions. The last three parameters 
are nuisance parameters accounting for foregrounds contributions to the SPT data set. 

lated to each other by a time-independent factor, as long as we assume that the annihilation 
cross section does not vary with temperature). A model-independent analysis, based on the 
expansion of Pann(^) in principal components, reveals that the CMB is mainly sensitive to 
the first principal component, peaking around z — 600, and at the next order to the second 
principal component, accounting for the redshift variation of f{z) around this same value. 

The goal of this section is to check these results with a simpler approach than a full 
principal component analysis. We will stick to the simple approximation for Pann(^) sug- 
gested above, involving two plateaus and one parabola. This family of functions has two 
free parameters, one amplitude and one curvature. We can choose to define the amplitude 
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Figure 5: Redshift dependent Pann(^), approximated as a two-parameter family of functions as 



described in the text, with 



,{z = 600) = 1- 10-^ m^s-^kg-i and a = -0.05, -0.1, -0.2. 



parameter pann at z = 600: 



Pann exp 
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2501 J 
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for z > 2500, 
for 30<z< 2500, 
for z < 30, 



(3.3) 



withpann = /(^ = 600) (ai;) /ttidm and o: < 0. With respect to the previous section, we now 
have a new dimensionless parameter a, that expresses the redshift dependence of Pann(^)- 
The question is whether this new parameter can be detected with current data: if not, the 
analysis of the previous section captures all the information that we can extract, with Pann 
standing for the value of the annihilation parameter near z ^ 600. 

We show in figure 6 the evolution of Xg and Tm for fixed pann and several values of a. 
As long as a remains small in absolute value {\a\ <^ 1), its impact is mainly on the slope of 
Xe{z) in the region z <^ z^ec- We expect this slope to be difficult to probe experimentally, 
since the CMB is mainly sensitive to the optical depth, which is an integrated quantity over 
redshift. 

When the redshift dependance of Pann(^) increases with a fixed normalization at z = 
600, the annihilation rate at high redshift increases. We expect to reach such large values 
that the decoupling time is not just slightly affected by DM annihilation, but radically 
postponed to a later time, because the massive energy injection from DM annihilation 
forbids hydrogen recombination. This happens for a < —3, as illustrated in figure 6. In 
this regime, the sound horizon at recombination is dramatically increased, and the CMB 
data will enforce a similar increase in the angular diameter distance to last scattering, in 
order to keep the same peak scale in multipole space. This will generate a correlation 
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between a and parameters such as the Hubble rate. However, the scale of the acoustic 
peaks and of Silk damping react differently to such a transformation, so we expect that a 
cannot be pushed to arbitrary negative values. This "extreme" regime could not be reached 
in the previous subsection: as long as we assumed a constant Pann, observational bounds on 
Pann prevented the annihilation rate to be too high around z ^ 1000. When comparing this 
model with CMB data, we first imposed no prior on a (apart from a < 0). We obtained a 
bound a > —5.3 (95% C.L.) and some non-trivial correlation between very negative values 
of a and other parameters. 




10° 10^ 10^ 10^ 10^ 10^ 10° 10^ 10^ 10^ 10^ 10^ 



1+z 1+z 



Figure 6: Free electron fraction and matter temperature Tm as a function of redshift for a 
constant Pann = 0, 1- 10~^ and 1- 10~^ m^s~^kg~^ (black and gray curves) and a redshift dependent 
Pann(^) with pann = 1' 10~^ m^s~^kg~^ and a = —0.1, —2 and 4, using recfast and assuming no 
reionization. 

However, this region in parameter space should not be taken seriously, because the 
realistic examples provided in [7] correspond to values of \a\ at most of the order of 0.1 or 
0.2. We performed a more "realistic" run with a top-hat prior —0.2 < a < 0. The results 
are summarized in the second column of Table 1. The data still gives no indication in favor 
of DM annihilation. The posterior probability of a is flat throughout the prior range, and 
the bounds on other parameters are essentially unchanged with respect to the run with a 
constant annihilation parameter, i.e. with a = 0. Even the two-dimensional probability 
contours in the space (^Pann, show no significant correlation between these parameters. 

These results are fully consistent with those of [10], showing that the first principal 
component peaks near z ^ 600. They also prove that current data is not sensitive to the 
second principal component, unless it has an unreasonably large amplitude like in the run 
with no prior on a. We could have defined our parameter pann at a different redshift: in 
that case, we would have expected to find a correlation between pann and a. The maximum 
of the first principal component can be seen as the "decorrelation redshift" between Pann 
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and a. 

In conclusion of this section, it appears that the moderate variation of f{z) (or equiv- 
alently Pann(^)) in the range 40 < z < 1000 suggested by the reahstic examples of [7] is far 
from being detectable with WMAP7+SPT data. Ref. [9, 11] reached a similar conclusion 
by comparing bounds on j^ann for some particular cases out of the possible f{z) functions 
presented in [7]. In the next section, it will be legitimate to neglect any variation of j>ann(^) 
at least until the redshift of halo formation. 



4. Annihilation in Dark Matter halos and Reionization 

Until now we considered that Dark Matter is uniformly distributed in the universe. It is 
well-known that structure formation generates a concentration of DM in halos that is likely 
to enhance DM annihilation. This issue has been widely discussed in the context of dark 
matter indirect detection in cosmic rays. It has also been pointed out that enhanced DM 
annihilation could be relevant for the reionization of the universe, and therefore for CMB 
physics [3, 4, 5, 8, 11, 12]. In this section, we wish to propose a systematic investigation 
of such effects, based on a generic parameterization of DM annihilation in halos, and a full 
parameter extraction from CMB data. 

4.1 Energy density release in DM halos 

The energy released by Dark Matter annihilation in halos can be written as (see e.g. [4, 12]) 

dE 

halos 



dVdt 



= P^(z)c^Pann, (4.1) 



where is the squared dark matter density averaged over space, that exceeds the square 
of the average dark matter density in presence of non-linear structures. In the halo model, 
this quantity is given by: 

dMj^i dr^^r^pl{r)\. (4.2) 

-^min \J ^ / 

Here Mmin is the minimal mass of DM halos, ^ the differential comoving number density of 
DM halos of mass M, r2oo the radius of a sphere enclosing a mean density equal to 200 times 
the background density, and ph the spherical DM halo density profile. The shape of density 
profiles is still a subject of controversy. If we consider for instance a Navarro- Frenk- White 
(NFW) profile [41], we can express the last integral as 



dr4nr p^{r) = {-q^ I /NFw(ch), (4.3) 

where z-p is the redshift of halo formation, p{zy) — 200/Ocf2M(l + ^f)^ the average mat- 
ter density within a radius r2005 and /new a function of the so-called halo concentration 
parameter c^- We recall that the critical density pc and density fraction parameters f2j 
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are defined today. In order to get an analytic approximation to dn/dM^ one can use the 
Press-Schechter formalism [42], leading to 

dn . pMd\iicr~^ . ,^ 

where pM = Pc^M is the average matter density today, and /(cr) the differential mass 
function. The variance of the linear density field cr(M, z) is given as usual by 

/•oo 

a2(M,z)= / P{k,z)W^{k,M)k^dk, (4.5) 

with P{k, z) the linear power spectrum at redshift and W{k^ M) the window function. 
Assuming that the collapse of the high density regions can be described by a spherical 
model, one can use a top-hat filter for W [43]. For the differential mass function, we could 
rely on the original function of Press and Schechter 

with 5sc — 1.28, or the more accurate function proposed by Seth and Thormen [44]. We 
could compute these terms exactly within the Boltzmann code, but the CMB spectra are 
not highly sensitive to the details of the halo model: they can only provide constraints on 
integrated quantities. Hence, it is irrelevant to search for high accuracy in this context. 
Instead, it would be very useful for the purpose of fitting CMB data to derive a simple, 
approximate parametric form for the energy injection function. To start with, we can use 
the fact that in a universe dominated by matter (i.e. any time between decoupling and 
z ^ 1), the redshift dependence of the variance a is somewhat trivial: 

a(M,z) = a(M,l)-^. (4.7) 



If we recollect all terms, we get the following contribution to the energy rate due to Dark 
Matter halos 



dE 



dVdt 



f^DM C^Pann(^) (1 + zf ^(1 + 2:F)^/NFw(Ch) 

halos 

- ^^^J dlnc7-HM,l) 2 ^.e(l + ^) ^^^p/ 51(1 + 



M„ 



dM v^2v^(t(M,1) V 8ct2(M,1 



(4.8) 



where we used the original Press-Schechter differential mass function for simplicity. The 
redshift-dependent integral simplifies with the change of variable u — r^ArTTT- If we 

2v 2cr(-M,l j 

define u^^z) = ^J^^^^^j^^^^y it reduces to 

/ (in— ^ exp (— n^) = erfc(nmm(-2)) , (4-9) 
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where erfc(x) is the complementary error function. It is suppressed at high redshift, before 
halo formation, i.e. as long as Umm{z) ^ 1. At low redshift, we do not expect the function 
Pann(^) to vary significantly, as can be seen in figure 4 of [7] for several examples. Hence, 
we can replace it by a nearly constant value Pann(O). In this case, the energy rate from 
annihilation in halos can be expressed as a function of only two parameters (beyond Pc^dm), 



dE 



dVdt 



p,2f^V'(l + ^f/herfc 



halos 



l + z 

1 + 



where is the characteristic redshift at which halos start to contribute^, 



2V^ 

Zh = —r 0-(M„ 



1-1 



(4.10) 



(4.11) 



and /h is a parameter related to the formation redshift and concentration of halos, and to 
the DM annihilation function at low redshift, 



fh = + ^F)^/NFw(Ch)j>ann(0) . 



(4.12) 



The parameter /h shares the same units as Pann- In principle, and the ratio fh/PanniO) 
should be inferred from a fit of the Press-Schechter formalism to detailed simulations of 
structure formation. However, there is no full consensus yet on the dynamics of halo 
formation and on halo density profiles. Moreover, these parameters should have a strong 
dependence on cosmological parameters, and also on the matter power spectrum at large 
A:, which is poorly constrained by observations. Hence we will treat Zh and /h as free 
parameters in our analysis. 

4.2 Beyond the on-the-spot approximation 

Sticking to the on-the-spot approximation, we could express the net ionization rate per 
interval of redshift + Ix^ using eqs. (2.7,2.8). The equation of evolution for Xg and Tm 
would receive extra contributions 



dxe 
dz 

dTyi 



dE 



dz 



{l-^ z)H{z)nH{z) dVdt 
1 dE 
' {1 + z)H{z)nH{z) dVdt 



1-Xe{z) f^^l-C{z) 



Er 



2 l + 2xe(z) 



1 



1 + fHe + Xe{z) 



(4.13) 
(4.14) 



with -^r^{z) given by the sum of the smooth density and halo density contributions. How- 
ever, this approximation becomes inaccurate at low redshift, as we shall see below. Well 
after recombination, the energy injection responsible for ionization and heating at a given 
redshift z comes from the decay products of DM annihilation at all redshifts z^ > z^ taking 
into account the fact that particle energies are redshifted between z^ and z, and that a 
fraction of the particles created at z^ interact along the line-of-sight and do not play a role 
at z. Moreover, in general, the interaction cross-section between z^ and z depends on the 



^In fact, the functions erfc(x) starts to raise at x < 2, so halos contribute below z < 2zh- 
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energy of each particle. Hence, the ionization rate obeys to a rather comphcated equation 
involving two integrals: one over dz' ^ and one over the energy of the particles created at z' ^ 
and interacting with the plasma at intermediate redshift. 

However, ref. [4] found that at low redshift, most of the ionization and heating is caused 
by photons produced by the inverse Compton scattering (ICS) of charged particles resulting 
from DM annihilation over CMB photons. It was shown by the authors of [8] that taking 
mainly this process into account leads to the simplified expression: 

(4.15) 



dz 



CUT 7(2;) 



{l + z)H{z) 



dz' 



1 + z 



dVdt 



{l + z')H{z') \l + z' ^ 

where e'^^^'^'-' is an absorption factor: it represents the fraction of photons produced around 
z' by ICS that already interacted with the inter-galactic medium and deposited their energy 
before z. Hence k can be approximated as 



k{z, z') 



car 



r' -dz" 

Jz {l + z")H 



(z") 



nniz") 



(4.16) 



not to be confused with the optical depth of CMB photons, featuring an extra factor Xe{z"). 
Finally, 7(2;) is a short-cut notation for 

-Xeiz) ( 1 , l-C{z) 



1 ^ 1 

3 \Ei Ea 



The integral in ^(z, z') can be performed analytically: 



k{z, z') 



2 nniO) 



Let us define a = cgt 



(l + ^)3/2_(i + ^/)3/2 ^ 

= . We notice that 

r 

dz'aVTT^e"^'^''^ 



-CUT 



_ H{z) H{z') _ 



1 



(4.17) 



(4.18) 



(4.19) 



Hence, the function A(z, z') = ay/\ + z' e'^^^'^'^ peaking in z' = z can be approximated 
with the Dirac function 6{z — z') in the limit in which it decreases with z' much faster than 
any other function in the integrand of equation (4.15). Writing (4.15) as 

3 



dz 



7(2) 



{l + z)H{z) 



dz' A(2,z') nniz') 



1 + z 
1 + z' 



dE 
dVdt 



(4.20) 



we see that in the approximation mentioned above, one recovers exactly the on-the-spot 
expression of eq. (4.13). But in the general case, we have to deal with the full integral. This 
is mathematically equivalent to keeping expression (4.13), with the on-the-spot energy rate 
replaced by an effective one, defined as 

3 



dE 



dVdt 



efr 



Jz nniz') V 



dE 



1 + z' J dVdt 



dE 
dVdt 



(4.21) 



- 19 - 



Similarly, this effective energy rate should be used in equation (4.14) to get the correct 
temperature evolution. 




1 10 100 1000 

1+z 



Figure 7: Effective energy injection rate rescaled by (1 + z)^, computed with or without the 
on-the-spot approximation, for either /h = 10000 or 1000. Other annihilation parameters are fixed 

to Pann = 10~^ and Zh = 20. 

We recall that the energy injection rate dE/dVdt includes a redshift dependent function 
f{z) accounting for the fraction of the energy released by DM annihilation that is actually 
absorbed by the gas, through ionization and heating processes. At low redshift, when 
annihilation is enhanced by halo formation and may contribute to the reionization of the 
universe, equation (4.21) refers specifically to ICS processes. Hence, for z ^ 100, the 
fraction f{z) just represents the fraction of the energy released by DM annihilation that 
goes into decay products with electromagnetic interactions (namely, photons and all charged 
particles), called /em in ref. [12]. 

4.3 Effects on Xg, Tm and the CMB spectra 

In figure 8, we compare the effect of DM annihilation in halos with that of the default 
reionization model implemented in CLASS and CAMB, based on a hyperbolic tangent centered 
in Zreio- The effect of DM annihilation on the ionization fraction is found to be very 
similar, except that it induces a slower reionization. The parameter controls the onset of 
reionization from halos, while /h controls its amplitude. For large enough values of /h, DM 
annihilation in halos can entirely reionize the universe before the current epoch, as shown 
previously in [4, 8, 12]. With the default reionization model, the ionization fraction Xg is 
larger than one at low redshift, because Helium reionization is also taken into account. In 
our model for DM annihilation in halos, we neglect helium reionization for simplicity, so 
that Xe is smaller or equal to one by definition. 

When including the effect of DM annihilation in halos, we work with the RECFAST 
mode of HYREC. Indeed, with RECFAST, we experienced numerical instability issues: the 
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Figure 8: Free electron fraction and matter temperature for Pann = 0, 10~^ and 10~^ m^s~^kg~^ 
(from bottom to top) and different values of /h and z^^ compared to the usual results for pann = 
and a single-step model for reionization from stars. All curves were obtained using hyrec in mode 
RECFAST. 



free electron fraction explodes and oscillates very rapidly already for small value of our 
parameters /h and z^- With HYREC in FULL modes, the only problem is that for large 
values of and /h, the ratio Tm/T^ may exceed one, falling outside the range of one 
interpolation tables. The RECFAST mode of HYREC is always well behaved. 

The right plot in figure 8 shows that the matter temperature increases a lot due DM 
annihilation in halos. Note also that for extreme values of the temperature Tm > 2 x 10^ K, 
using RECFAST's case-B recombination coefficient becomes inaccurate [45]. We will see 
anyway in section 4.6 that such large values are in contradiction with constraints on the 
temperature of the inter-galactic medium at z < 4, as inferred from Lyman-a observations: 
this will provide an addition constraint on the DM annihilation rate. 

The signature of DM annihilation on the primary CMB anisotropy spectrum is found 
to be very similar to that of reionization. In addition to the peak shifting and damping 
due to a non-zero pann parameter, the halo effect controlled mainly by /h leads to an 
overall suppression of temperature/polarization power for / > 30, and an enhancement of 
polarization for / < 30. We can anticipate that the CMB alone can hardly discriminate 
between the contribution of reionization from stars and from halos, since the CMB spectra 
probe mainly the optical depth, i.e. the integral of Xe over time. However, the fact that 
DM induces a slow reionization process starting at high redshift^ implies that the step- 
like suppression of temperature and the low-/ polarization bump are smoother and wider 
than with the default reionization model. To illustrate this, we compare in figure 9 the 
low-/ polarization spectrum for two models with the same optical depth. Accurate CMB 
polarization data limited only by cosmic variance on large angular scale may probe such a 

^In the CMB analysis of the next subsections, Zh is found in the range from 20 to 30, implying that 
halos start contributing between 40 and 60, well before star formation. 
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difference. 



Standard reionization 
Reionization by DM halos 
Reionization by DIVI halos and stars 




Figure 9: Low-/ polarization spectrura for the three best-fitting models assuming reionization from 
stars (with the usual single step parameterization), from DM annihilation in halos, or from both 
with an additional Gunn-Peterson prior. 

4.4 Can Dark Matter annihilation alone explain reionization? 

We wish to check whether WMAP7 and SPT data are compatible with the assumption that 
the reionization of the universe can be explained entirely by DM annihilation in halos, as 
suggested in [4, 8, 12]. It is rather obvious that the free parameters and /h of our model 
can be adjusted in such way that the reionization optical depth is be compatible with the 
WMAP7 best-fitting value. However, we have seen that DM annihilation can only induce 
slow reionization starting at high redshift, and induce a wider step (resp. bump) in the 
low-/ temperature (resp. polarization) spectrum. A priori, this may lead to a value of the 
maximum likelihood significantly lower for the annihilation model. In addition, an analysis 
with free and /h could lead to preferred values in strong contradiction with expectations 
from structure formation and halo models. 

The results of our MONTEPYTHON run with WMAP and SPT data are summarized in 
the third column of Table 1 and in the triangle plot of figure 10. The new free parameters 
and /h are not degenerate with other parameters, so the credible interval for the usual 
ACDM parameters and A^eff are unchanged with respect to the standard model without 
annihilation. There is instead a significant correlation between and /h: if halos form 
very late, a very large amplitude parameter /h is needed in order to get the same optical 
depth. The effective chi square Xgff — ~2 InC is higher for DM reionization than single-step 
star reionization, but only by 0.8, showing that the data shows no strong preference for one 
model against the other. 

The characteristic redshift Zh is found in the range 11 < < 40 (95%C. L.). This 
parameter has a strongly non-gaussian posterior probability, with a mean value of 23, but 
a best-fit value of 29. The shape of the erfc(x) function is such that the halo contribution 
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starts raising around z ^ 2z\^. Values of close to 20 imply a raise in the range 40-50, 
which is plausible from the point of view of structure formation. 

To see whether the required value of halo concentration is sensible, we need to make an 
assumption about the DM annihilation amplitude, since the data is sensitive to /h, i.e. to 
the product /NFw(ch) x J>ann(0). If we first assume a value of pann at z ^ 600 saturating our 
CMB bound, Pann(600) ^ 4 x 10"'', we expect that at low redshift this parameter will fall 
to approximately j>ann(0) ^ 10~^. Then the best-fitting value /h ^ 12600 m^/s/kg requires 

/NFw(ch) - _^^^^3 ^ - 8200 , (4.22) 

200(1 + zyYv^^A^) 

where we also assumed ^ 60. The quantity /NFw(ch) is poorly constrained, but models 
of halo formation suggest an order of magnitude ranging from 10^ to 10^ Hence the " reion- 
ization from DM annihilation" model points towards a reasonable value of the concentration 
parameter. If /NFw(ch) is of the order of 10"^, then constraints on DM annihilation from 
the smooth background and from halos are comparable. If /NFw(ch) is of the order of 10^ 
(resp. 10^), then constraints from annihilation in the smooth background (resp. in halos) 
are stronger. Indeed, the bound on Pann coming from annihilation from halos can be found 
using the relation 

P.„„(. = 600) = 3.3 X 10-" "\3, , A . (4.23) 

5 (1 + ^f)'^ /NFW^(Ch) 

Taking /^^^ = 24500 m^/s/kg (95% C.L.) from our analysis, this implies 



Pann(^ = 600) < 0.81 X 10" 



-6 



j9ann(600) 61^ 10^ 
. 5pann(0) (1 + Z^f /NFw(Ch)_ 



mVs/kg (95%C.L.) (4.24) 



The low-/ polarization spectrum for the best fitting model is shown on figure 9. Given 
its distinct shape due to an early and slow reionization process (with respect to star reion- 
ization), we expect future small-scale polarized measurements by Planck and other CMB 
experiments to improve the bound on /h. 



4.5 The Gunn-Peterson effect 

It was realized in 1965 by Gunn and Peterson [46] that the observation of redshifted Lyman- 
a absorption lines in quasar spectra was a very sensitive probe of the presence of neutral 
hydrogen along the line of sight, and hence of the ionization fraction of the universe at 
different redshifts. Since even a small fraction of neutral hydrogen leads to a clear signature, 
we have some evidence that the universe was almost fully ionized until z ^ 6, since quasars 
at such a redshift show a very small level of Lyman-a absorption. More precisely, according 
to [47], the fraction of neutral hydrogen xhi has to satisfy: 

• for z > 6, xhi > 10~^ (xhi might even be equal to 10"-*^), 

• for z < 5.5, xhi < 10~^, 
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Figure 10: Assuming a ACDM model with free electron fraction, dark matter annihilation (back- 
ground and halos) and no extra reionization from stars, marginalized probability distribution of the 
annihilation parameters Pann, /h and given WMAP7 and SPT data. 

(see also [48] for a recent constraint at z 2^ 7). Thus there seems to be an abrupt transition 
between z = 5.5 and z = 6. This raises some tension with the simplest model of single- step 
reionization from stars, in which the ionization fraction is assumed to evolve roughly like a 
hyperbolic tangent centered on a reionization redshift Zreio- The problem is not related so 
much to the precise shape of the step, but to the fact that any abrupt step should be centered 
near z = 6 or 7 to comply with Gunn-Peterson observations, instead of Zreio — 10.6 =b 2.4 
(95%C.L.) to explain the optical depth Treio = 0.088 zb 0.015 measured by WMAP [37]^ A 
single-step reionization at z ^ 10 or even 8 would lead to xhi < 10~^ at z = 6. 

The model of the previous subsection, in which reionization is caused entirely by DM 
annihilation in halos, also fails to explain Gunn-Peterson observations for the opposite 
reason: reionization is then so slow that all allowed models have xhi > 10~^ at z = 5.5. 

There could be several solutions to this problem: 

• the Gunn-Peterson bounds may be wrong or not correctly interpreted. These bounds 
are in fact model-dependent and controversial, since they rely on assumptions con- 
cerning the density and temperature of the inter-galactic medium, and the ultra-violet 
background. Observations at z ^ 6 could be explained with alternative models for 
the IGM and UV background, instead of incomplete reionization [49]. The evidence 
that the universe is fully ionized below z ^ 5.5 could also disappear with different 
assumptions, for instance in the context of inhomogeneous reionization [50]. 

^To be more precise, the CMB constraint is dominated by the measurement of low-/ E-type polarization 
by WMAP, and depends on the assumed cosmological model 



- 24 - 



• the cosmological model describing our universe may have extra ingredients (not nec- 
essarily related to reionization) such that a good fit to WMAP data can be obtained 
with single-step reionization at Zreio ^ 6 or 7. 

• reionization may be caused by different population of stars forming at different red- 
shifts. The single-step model is too naive and should be replaced by a model with 
at least two steps. The late one should take place around z ^ 6 or 7 to account 
for Gunn-Peterson observations. The early one, possibly related to the generation of 
massive, metal-free stars [51, 52], should partially reionize the universe and enhance 
the optical depth. 

• reionization may be caused both by star and by the decay or annihilation of some 
particles. The possibility of enhancing reionization with sterile neutrino decay has 
been proposed by [53]. In the case of annihilating DM, ref. [8] suggested that DM 
annihilation in halos may start to slowly reionize the universe. At a redshift close 
to six, star formation processes take over and quickly ionize the remaining hydrogen 
atoms. 

In this subsection, we wish to test the last paradigm. It is a priori not obvious that any 
model of this type can work, because in order to explain the observed optical depth, DM 
annihilation may need to be so large that in any case xhi < 10~^ at z = 6. Fortunately, we 
will see that this mixed model nicely complies with Gunn-Peterson and CMB constraints, 
and points to plausible halo parameter values. 

We added a Gunn-Peterson prior to WMAP7 and SPT data and ran MONTEPYTHON 
again. More precisely, we impose two top-hat priors 10~^ < xhi(6) < 1 and < xhi(5.5) < 
10~^. We neglect Helium reionization for simplicity. In this approximation, Xe just repre- 
sents the fraction of ionized hydrogen, and we have xhi = ^—Xe- Our Boltzmann code CLASS 
simulates mixed reionization in the following way. For each model, the ionization fraction is 
first computed down to z = neglecting reionization from stars, using HYREC (in RECFAST 
mode) in order to avoid numerical instability. The effect of stars is then implemented "by 
hand": below some arbitrary redshift Zreio, Xe{z) is cut and matched continuously to a 
half- hyperbolic tangent centered on Zreio, reaching an asymptotic value of one for z ^ 
(see one example of such models in figure 11). The precise shape of this function is in fact 
identical to that for ordinary single-step reionization in CLASS and CAMB, except that only 
the side z < Zreio of the step-like function is used, and that the transition width parameter 
is decreased to 6z = 0.2 in order to model a very fast reionization process. 

Our results are summarized in the last column of Table 1 and in the triangle plot of 
figure 12. With respect to the previous model of section 4.4, we have one more parameter 
Zreio, that is very constrained by the Gunn-Peterson prior: it can only fluctuate in the range 
6.4 < Zreio < 6.8 (95% C.L.). The posterior probability of /h and are shifted to slightly 
smaller values, since DM annihilation in halos is only expected to contribute to a fraction 
of the optical depth. For all other parameters, the results are essentially identical to those 
of the previous case. The minimum effective chi square is also unchanged. The discussion 
of section 4.4 concerning /h and still applies: is fixed to a range that coincides with 
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Standard reionization 
Reionization by DM halos 
Reionization by DIVI halos and stars 




(Uz) 

Figure 11: Free electron fraction for the three best-fitting models assuming reionization from stars 
(with the usual single step parameterization), from DM annihilation in halos, or from both with an 
additional Gunn-Peterson prior. All curves were obtained using HYREC in mode RECFAST. 

expectations from structure formation, and values of /NFw(ch) in the range from 10^ to 10^ 
could be accommodated provided that pann and /NFw(ch) fulfill the relation (4.23) with /h 
in the range 2100 < /h < 28400 (95% C.L.). 

In conclusion of this subsection, we see that this mixed model for reionization is interest- 
ing: DM annihilation in halos could explain the value of the optical depth probed by CMB 
data, while reionization from star formation at z 2^ 6.5 would complete the reionization 
process and explain Gunn-Peterson observations. 

4.6 Including an upper bound on the IGM temperature 

The best-fitting models of sections 4.4 and 4.5 have a halo parameter of the order of 
lO^m^/s/kg, leading to a matter temperature of the order of lO^K at low redshift (see 
figure 8). This estimate of the average matter temperature in the universe should be taken 
with a grain of salt, since we did not account for inhomogeneities in the matter distribution, 
nor for the thermodynamical evolution of the inter-galactic medium (IGM) during structure 
and star formation. 

However, Lyman-a observations suggest that the IGM temperature is of the order 
of a few times lO^K in the redshift range 2 < z < 4.5. Reference [5] pointed out that 
these measurements should provide at least an upper bound on the average temperature 
enhancement due to DM annihilation. 

In other words, the results of the previous two subsections are compatible with Lyman- 
a observations only if we are modelling the temperature evolution incorrectly. The ansatz 
that a fraction (1 + 2xe)/3 of the energy injected into the gas by DM annihilation goes 
into heating might be incorrect at low redshift; or the IGM temperature growth might be 
limited by some temperature regulation mechanisms not described by our simplistic set 
of equations. If instead our temperature evolution law is realistic, then DM annihilation 
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Figure 12: Assuming a ACDM model with free effective neutrino number, dark matter anni- 
hilation (background and halos) and single-step reionization from stars, marginalized probability 
distribution of the annihilation parameters Pann, /h, and Zreio given WMAP7 and SPT data and 
an additional Gunn-Peterson prior. 

cannot explain the reionization of the universe alone, and cannot even contribute sufficiently 
to reionization at z ^ 6 in order to explain Gunn-Peterson bounds with a mixed reionization 
model, based on annihilation plus a single-step star formation process. 

It is still interesting to perform a parameter extraction with a prior on Tm, while assum- 
ing a mixed reionization from halos and stars, in order to check whether IGM temperature 
estimates provide a stronger bound on pann than the CMB alone. In this section, we will re- 
peat our analysis with a conservative upper bound on the IGM temperature at low redshift, 
inspired from figure 6 (left) in [54]: 

Tm{z = 2) < 3.2 X lO^K. (4.25) 

We implemented this constraint in the form of a top-hat prior in MONTEPYTHON. We 
checked in presence of such a prior, the model of section 4.4 in which reionization is caused 
entirely by DM annihilation requires unrealistically high values of z/^ ^ 100, totally incom- 
patible with structure formation models. 
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Figure 13: Assuming a ACDM model with free affective neutrino number, dark matter anni- 
hilation (background and halos) and single-step reionization from stars, marginalized probability 
distribution of the annihilation parameters Pann, /h, and Zreio given WMAP7 and SPT data, a 
prior 20 < < 30 and an upper bound on the IGM temperature. 



Like in the previous subsection, we implemented star reionization into CLASS "by 
hand": below some arbitrary redshift Zreio, Xe{z) is cut and matched continuously to a 
half- hyperbolic tangent centered on Zreio, reaching an asymptotic value of one for z ^ 0. 
We kept the transition width parameter at its default value: 5z — 1.5. We also imposed 
a top-hat prior 20 < < 30 in order to ensure that halos form at a realistic redshift, 
compatible with simulations of structure formation. 

Our results are summarized in the last column of Table 1 and in the triangle plot of 
figure 13. The upper bound on /h is reduced by one order of magnitude due to the IGM 
temperature constrain: the 95%CL upper limit on /h decreases from /^^^ = 24500 m^/s/kg 
to 1200m^/s/kg. This leads to a stronger bound on Pann- 



p^^^{z = 600) < 0.04 X 10" 



Pann(600) 61^ 



10^ 



mVs/kg (95%C.L.) (4.26) 



_ 5;?ann(0) (1 + ^f)^ /NFw(Ch) 

Assuming that the factor between brackets is equal to one, this bound is ten times stronger 
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than the one inferred from annihilation in the smooth component only. 

The heating effect of DM annihilation may also enhance the kinetic Sunyaev-Zel'dovich 
effect [55], leave a signature in secondary CMB anisotropics, and provide a further test of 
this model [12]. We do not study this aspect in our work. 

5. Conclusion and Outlook 

We studied different possible contributions of annihilating Dark Matter to the thermal his- 
tory of the universe. We confirmed previous results that the annihilation of the background 
DM distribution has non-trivial effects on the CMB, leading to the constraint 

P^nn{z - 600) < 0.46 X 10"VVs/kg (WMAP7 + SPT, 95%C.L), (5.1) 

with a negligible impact of the variations of Pann(^) in the range 100 < z < 2500 suggested 
by a realistic study of DM annihilation channels. 

We also showed that DM annihilation in halos could explain entirely the reionization 
of the universe from the point of view of CMB observations. In addition, if the constraints 
^Hi(6) > 10~^, xhi(5.5) < 10~^ inferred from the Gunn-Peterson effect hold, and if we as- 
sume that reionization from stars takes place abruptly in one step, then a mixed model with 
DM annihilation in halos and star formation at z ^ 6.5 could explain simultaneously CMB 
observations and the above bounds. However, these models tend to reheat the IGM well 
above the typical temperatures indicated by Lyman-a observations, unless our modeling of 
the matter temperature evolution at low redshift is incorrect. 

Our most important conclusion is that constraints on DM annihilation in halos tend 
to be stronger than those from the smooth background distribution of DM, especially if 
we include a realistic upper bound on the matter temperature at low redhsift. Assuming 
Press-Schechter theory and NEW profiles, we see from eq. (4.26) that for /NFw(ch) — 10^, 
zy = 60 and [j>ann(600)/j)ann(0)] = 5, the constraint coming from halos and from the smooth 
background are roughly equal. If in reality halos are more concentrated than in this simple 
model, then constraints on Pann from annihilation in halos superseed those from annihilation 
in the background. We summarize our constraints on j>ann and their implications for the 
DM mass and cross-section in figure 14. A WIMP with standard thermal cross-section 
(av) 2^ 3 X 10~^^ cm^/s is constrained by annihilations in the smooth background to have a 
mass larger than 35 x f(z = 600) GeV/c^ (95% C. L.). If halos are very concentrated and 
can be modeled by /NFw(ch) ^ 10^ (resp. 10^), the bound increases to about 200 GeV/c^ 
(resp. 2000 GeV/c^). 

Many aspects of this analysis could be improved. For example, we neglected Helium 
throughout the whole discussion. The effect of Hehum has been studied by Galh et al. [9] 
(neglecting halo effects) and does not change the result significantly. The energy fraction 
going into ionization Xi Lyman-a excitation Xa was also approximated, motivated by 
a common sense argument by Chen and Kamionkowski [32]. However, the exact behavior 
of these quantities have a negligible effect on the CMB. The average DM density during 
non-linear structure formation has been approximated with a basic Press-Schechter model 



-29- 




Figure 14: Summary of our bounds on pann translated into constraints for the DM annihilation 
thermal cross-section {(jv) and mass moM- Points in the shaded regions are above the 95% preferred 
region for Pann(^ = 600), considering only annihilation in the smooth DM background, and assuming 
either f{z = 600) = 0.2 (dark shade) or f{z = 600) = 0.9 (light shade): these two extreme 
assumptions cover the plausible range for f{z = 600) in the case of WIMP annihilation, see [7]. The 
three black lines correspond to the CMB bounds inferred from DM annihilation in halos, assuming 
/NFw(ch) = 10^ (top), 10^ (middle) or 10^ (bottom), and taking in all three cases = 60 and 
f{z ~ 0) =0.1 (or in other words, f{z = 600) = 0.5 and [pann(600)/pann(0)] = 5). When a realistic 
upper bound on the matter temperature at low redshift is taken into account, the bounds move to 
the green lines. The horizontal lines shows the standard WIMP thermal cross-section. 

and NFW profiles. We could have imposed priors on the parameters of this model in- 
ferred from N-body simulations, or tried different profiles (Einasto profile, etc.), or a more 
realistic differential mass function [44]. Instead of the Press-Schechter model, we could 
have accounted for halo formation using the excursion set formalism [56]. One could try 
to model the matter temperature evolution at low redshift more accurately, taking into 
account matter inhomogeneities and the thermodynamical evolution of the IGM. However, 
all these refinements are probably unnecessary at the moment, given the large error bars 
on the optical depth inferred from CMB observations. 

Throughout this work, we assumed that DM annihilates. A similar study can be 
performed in the case of decaying dark matter [13, 15, 16, 17]. In that case, the energy 
injection rate varies like pdm (instead of p^^)' ^^^^ + ^)^- Hence, the effect of DM 
decay in the smooth DM background is not very different from the effect of DM annihilation 
in halos, studied in section 4. Note however that for a wide range of masses, constraints on 
the DM lifetime inferred from current CMB observations are not as strong as those inferred 
from cosmic rays [57, 58]. 

In a few months from now, results from the Planck satellite data may lead to a signifi- 
cant improvement of these bounds, and bring complementary information on the DM mass 
and cross-section (or lifetime) with respect to direct and other indirect DM search. 
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Appendix A: modifications in HYREC 



In HYREC, the evolution equation are written in function of time for Xg and of In a for Tm- 
In addition, the units are CGS+eV for temperatures, except in the two functions describing 
the temperature evolution, where we have Kelvin. So, the equations are 



H 



dXe 



1 



din a H 



_E2i 



3nH 



— - - Zlu + 77 

a ma H 

2 l + 2Xe 

+ 



3meC 1 + Xe + fne 

1 



1^1- 

Ei E, 



c 



Tm) 

p2c2f22^(l + z)Vann(^) 



(5.2) 



(5.3) 



3^6 3nH 1 + Xe + fHe 

These equations are equivalent with the ones in RECFAST, but the C-factor is now defined 



as: 



C 



Pb + \RLya + \^ls,2s ' 



(5.4) 



with the two-photon rate A25,i5 = 8.22458 s ^ and the escape rate of Lyman-a photons 



Rl 



SttH 



Quasi steady-state equation 

In HYREC, we can find a function describing the temperature evolution in the quasi steady- 
state approximation. In general, as seen above, the equation for the temperature in the 
presence of annihilating Dark Matter is 



din a 



-2Tm + 7{Tr - Tm) + 



din a 



(5.5) 



DM 



The quasi steady-state approximation consists of considering the case when the second term 
in equation (5.5) is bigger than the other two, i.e. when 7^1. In this situation, Tm ^ Tr 
-Tm, thus 



and 1^ 

a In a 



Tm 



1 + 7 



Tr , _i dTu 
1+7 



din a 



(5.6) 



DM 
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The different modes of HYREC 



The above evolution equation for Xg is used when HYREC runs in the modes peebles or 
reef ast. For the EMLA mode, there exist two different C-factors, namely C2s and C2p 

C2s = (5.7) 

^2s - ^2s^2p 

A2sls 



RLya + R2p^2s 
-L 2p — ^2p^2s 



The inverse life times are defined by: 



= B2s + R2s^2p + A25,ls, (5.9) 
= + R2p^2s + ^Lya, (5.10) 

where Bi are the ionization coefficient and Ri^j the transition coefficients. We can take 
R2s^2p — 3i?2p^2s since there are 3 times more states in 2p than in 2s. 
C2s {C2p) represents the probability that a hydrogen atom initially in the 2s (2p) state 
reaches the ground sate before being ionized. The Lyman-a line is the excitation from Is 
to 2p. So the new factor in equation (5.2) should be^ C = C2p. Exactly the same approach 
is used for the full mode of HYREC. 
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